Three-dimensional theory of quantum memories based on A-type atomic ensembles 
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We develop a three-dimensional theory for quantum memories based on light storage in ensembles 
of A-type atoms, where two long-lived atomic ground states are employed. We consider light storage 
in an ensemble of finite spatial extent and we show that within the paraxial approximation the Fresnel 
number of the atomic ensemble and the optical depth are the only important physical parameters 
determining the quality of the quantum memory. We analyze the influence of these parameters on 
the storage of light followed by either forward or backward read-out from the quantum memory. 
We show that for small Fresnel numbers, the forward memory provides higher efficiencies, whereas 
for large Fresnel numbers, the backward memory is advantageous. The optimal light modes to 
store in the memory are presented together with the corresponding spin-waves and outcoming light 
modes. We show that for high optical depths such A-type atomic ensembles allow for highly efficient 
backward and forward memories even for small Fresnel numbers F 
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I. INTRODUCTION 

The natural way to transmit quantum information is to 
use photons. A photonic state, however, has to be stored 
locally in order to process the information or to use it 
at a later time. To this end, one needs an efficient and 
controllable quantum interface between light and matter 
that will store incoming light as a stationary excitation 
and release it at a later time preserving quantum correla- 
tions. This can be realized by, for instance, an ensemble 
of atoms or impurities in a solid state host, see Ref. [l| for 
an extensive review. The resulting quantum memory is 
an essential part of many quantum information devicej. 
It can allow for long distance quantum communication [21 
and, in connection with generation of entanglement and 
its swapping it might be a basic building block for a 
quantum internet 

One way to achieve strong coupling between light and 
matter is to use a single atom in a cavity [5|-[3|. This 
is, however, hard to realize experimentally and allows 
only for storing a single mode. An alternative approach 
is to use atomic ensembles, where the coupling strength 
scales with the number of atoms Its scalability al- 
lows for multimode memories @, [s] and it is easier to 
implement in practice. Quantum memories have already 
been demonstrated in a number of experiments based on 
atomic ensembles, see e.g. as well as solid state 

systems, see e.g. [l7l - [l9| . Despite the fast experimental 
progress in the field of single and multimode quantum 
memories, little work has been done on describing the full 
spatial profile of the excitations d, [20l - [23 [ . The widely 
used one-dimensional theory for spatial quantum mem- 
ories assumes only a single transverse mode [l^, [l^l . It 



has been shown that this theory is correct for infinitely 
large atomic ensembles [l|, H, . 

Within the one-dimensional approach [2^, [23| it has 
been demonstrated that the optimal control strategy for 
storage and retrieval of light can be applied to a wide 
class of quantum memory schemes including electromag- 
netically induced transparency (EIT'l [25| . off-resonant 
Raman (2^, and photon-echo based t27| methods. All 
of these schemes yield identical maximum efficiencies. 
Moreover, the maximal efficiency depends only on one 
physical parameter: the optical depth of the medium d. 
The retrieval of the light stored in the atomic spin- wave 
can be interpreted in terms of constructive interference 
of light emission in the forward direction set by the di- 
rection of the collinear control beam giving an effective 
rate 7^. Thus the efficiency is given by the ratio between 
this emission and the unwanted one: 
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The efficiency is independent of the detuning and the 
temporal shape of the control field and only depends on 
the shape of the stored quantum light pulse and the op- 
tical depth d. Furthermore, the ideal storage process can 
be seen as the time reversed of the read-out, thus also 
depending only on the optical depth d. 

In this paper, we develop a three-dimensional theory 
of quantum memory for light based on A-type atomic en- 
sembles, where two long-lived atomic ground states are 
used. We solve the problem of light storage in an atomic 
ensemble of finite spatial extent, where the spatial distri- 
bution of the stored stationary excitation (the spin- wave) 
depends on the transverse profile and temporal shape of 
the input light pulse. It is shown that within paraxial ap- 
proximation there are only two physical parameters that 
determine the quality of the memory: the optical depth 
d and the Fresnel number of the atomic ensemble F, 
whereas in the one-dimensional theory, where F ^ 00, 
only the optical depth matters. The same conclusion 



was reached in Ref. j20|, where a three-dimensional the- 
ory of the closely related problem of stimulated Raman 
scattering was developed. There, semi-analytical expres- 
sions for total emitted light were derived, whereas in this 
paper we perform exact numerical calculations describ- 
ing the full process of storage followed by retrieval of 
light. In addition, we optimize this combined process 
to obtain the best performance of the spatial quantum 
memory and its dependence on the physical parameters 
F and d. We calculate the efficiencies of storage followed 
by either forward or backward retrieval of light and study 
their dependence on the two crucial parameters of the 
atomic ensemble. We find that the direction which has 
the highest read-out efficiency depends on the Fresnel 
number of the ensemble. The optimal spin-wave modes 
with highest achievable efficiencies of the combined pro- 
cess of storage followed by retrieval from the memory are 
presented together with the optimal incoming and out- 
going light modes. In addition, we show that the optimal 
input light pulses are well approximated by the product 
of a temporal and transverse spatial shape. This is ad- 
vantageous from an experimental point of view, since one 
does not need to vary the transverse profile in time. We 
show that for time independent driving the optimal out- 
put light modes are the complex conjugates of the time 
reversed input modes. In addition, we calculate the rele- 
vant experimental parameters of the optimal input light 
modes, that is the beam waist and the distance between 
the focal plane and the beginning of the ensemble. As for 
the one-dimensional large class of resonant (EIT) 

and off-resonant Raman based experiments are contained 
in the presented three-dimensional theory. We show that 
for optically dense atomic ensembles, the studied quan- 
tum memory leads to high efficiencies of storage followed 
by backward or forward read-out even for small Fresnel 
numbers of the sample {F > 0.1). 

The remainder of the paper is organized as follows: In 
Sec. ini the model of the quantum memory is presented. 
Next, in Sec. IIIIl we analyze the read-out process and 
its dependence on the crucial physical parameters: the 
optical depth and the Fresnel number of the ensemble. 
The adiabatic storage and retrieval is studied in Sec IIVI 
where the analytical expressions are derived in Sec. lIV Al 
The results of numerical calculations for both read-out 
directions are presented in Sec. IIVBI together with the 
discussion of the implications of the three-dimensional 
theory. Sec. |V] concludes the paper with final remarks. 
In addition, some technical details and further analysis 
are presented in the Appendixes. 



II. QUANTUM MEMORY MODEL 

We consider an atomic ensemble of finite spatial extent 
which contains atoms with a A-type level structure [see 
Fig. 1 (left)]. The weak quantum field carries the infor- 
mation to be stored and couples states |0) and |e) with 
coupling strength g. The strong classical field couples the 




FIG. 1. (Color online) (left) A schematic plot of the level 
scheme of the considered A-type atoms with two long-lived 
ground states |0) and The quantum field carries the in- 
formation to be stored and couples |0) and |e) with coupling 
strength g. States |1) and |e) are coupled by the strong clas- 
sical control field with Rabi frequency Q{t). A denotes the 
detuning of the light from the excited state |e), which can 
spontaneously decay at a rate 7. (right) Atomic ensemble 
of cylindrical symmetry with Gaussian distribution of atoms 
in the radial direction. The density along the z direction is 
constant. The length of the sample is L and the width is 
characterized by a±. 



states |1) and |e) with Rabi frequency n{t) and controls 
the storage and retrieval of light into and from the mem- 
ory. The storage is performed as follows: from each exci- 
tation removed from the quantum light field one atomic 
excitation goes from |0) to |1) via the excited state |e) 
and the information is stored in a collective state of the 
ensemble spin- wave excitations S ~ |0)i(l|. 

We assume an atomic ensemble with cylindrically sym- 
metric density distribution [see Fig. [T] (right)] and the so- 
lution strategies presented below will work for any such 
atomic density distribution. For concreteness, however, 
we choose a Gaussian distribution in the radial direc- 
tion with n{r) = no exp ^~2a^)' ^here the normaliza- 
tion factor is no = Nj^/{2'nLG']_) with A'a — J (P'rn{'r) 
the number of atoms, and a±_ the width of the ensemble 
(cr^ being the corresponding cross section area) . For sim- 
plicity, we assume a constant density along the z axis, no 
for z G [0,L], so that our ensemble forms a cylinder of 
length L. 

To solve the three-dimensional problem of storage of 
a set of light modes, we expand the electric field into a 
complete and orthogonal set of modes fk.mn(r), which are 
solutions to the Helmholtz equation. For simplicity, we 
employ the paraxial approximation here and assume that 
the polarization Ck factors out fk,mn{r} = ekfk.,mn{r). 
The light beams under consideration are cylindrically 
symmetric, collinear, and propagate in the z direction 
and are conveniently described in a basis set of Bessel 
beams: 



fk.mniP, 4*1 ^) — NmnS^"^'^ Jmik_L,mnP)t 



(2) 
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where we have factored out the axial position z. Here, 
^mn = l/yTri?^ J^_|_i(Amn) is the normahzation factor, 



n-th zero of the Bessel function of order m of the first 
kind Jm, ko = ojo/c is the wave vector of the quantum 
hght with central frequency loq and speed c. Since we are 
working in the paraxial regime, where the perpendicular 
wave numbers are small, k±^rnn ^ k, we make the ap- 
proximation fc||,m„ - fc w -k±^rnn/{2k) « - A;_L,m„/(2A;o). 
Here, we have also imposed the boundary condition 
fk,mn{p = R,4>iz) = that all fields vanish at a vir- 
tual cut-off radius R. Thus we obtain a discrete set of 
transverse mode functions Umn{r±) that is complete and 
orthogonal 



= Jk"^ - k^, , — , Xmn is the 



(3) 
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In principle, one can use an arbitrary basis, but we choose 
this particular basis where one can factorize the depen- 
dence on the longitudinal z and transverse position f±, 

fk,ran(r) Wm„(rl)exp i {k- z . As we will 

show below, this makes the numerical calculations of the 
problem manageable. 

The expression for the; c;k;(;tric field of the quantum 
light expanded on this basis set fk mn{r) is then given 

by 



mn,k 



where I is the length of the quantization volume and eq 

is the polarization of the quantum field. We assume that 
both quantum and classical light fields arc narrowband 
fields centered at coo = Weo — A and = ^ei — A, respec- 
tively, so that iOk « wq. Here, A is the detuning from the 
excited state for both light fields with the corresponding 
atomic transition frequencies Weo and ujei ■ We introduce 
slowly varying light operators for the quantum light field, 
which depend on time and position 



S*'^°*. (5) 



The annihilation operators of the light field amn{z,t) are 
the expansion coefficients into the transverse mode basis 
and have the commutation relation 

Qmn (z) , al^,^, {z') = cSmm' ^nn' 6{z - z'). (6) 

With the chosen normalization, where c appears in the 
commutator, the photon intensity in mode mn at the 
position z is given by al^j^{z)amn{z)- The electric field 
rewritten by means of these new light operators. 
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+ H.c] 



is now an expansion on the set of transverse mode func- 
tions Um„(rj_). 

We also introduce slowly varying operators for the 
atomic quantities of interest, which depend on position 
and time: 

1 =<5(=')(r-rl)aS (7) 
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P{f,t) 
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^(3)(f_f.)CT«e-*'"°*+''=°^ (9) 
5^'\f-n)a'i}. (10) 



Here, k^ — Wq/c is the wave number of the classical con- 
trol field and a^^, = are the internal state op- 
erators of the i-th atom. S{f,t) denotes the collective 
spin-wave excitation in the atomic ensemble, P{f, t) is 
the polarization/coherence between the ground state |0) 
and the excited state |e), <Tei(r, describes the coher- 
ence between the second ground state |1) and the excited 
state |e), and aee(r, i) is the occupation of the excited 
state. For weak excitations below saturation, this ex- 
cited state operator can be neglected. We also assume 
that the atoms are stationary. The normalization of the 
atomic operators is chosen so as to achieve the same-time 
commutation relations are 



Sir,t),S\r',t) 
P{r,t),S^{r',t) 



P{f,t),P\r',t) =6^^\f-f{ll) 
1 ^^2) 



n{r) 



It will be convenient to express the atomic operators in 
the same basis as the light operators. Therefore, we also 
expand the operators for the spin- wave and the polariza- 
tion on the transverse mode basis 



S{f,t) 
P{f,t) 



E 

mn 



^mni^ -i-} ^"in{z , t) , 
{f±)Pmn{z,t). 



(13) 
(14) 



Having defined and expanded the slowly varying 

atomic and light operators, we now consider the Hamil- 
tonian for the system. The free light Hamiltonian reads 



E 

k.7nn 



(15) 



with the linear dispersion relation cok = ck. The atoms 
in the ensemble are described by the Hamiltonian 



(16) 



where the transformation to the rotating frame has been 
performed with respect to the laser frequency luq- 
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In order to obtain the Hamiltonian describing the in- 
teraction of hght with the entire atomic ensemble, we 
start out by considering the coupling of light and a sin- 
gle atom in the dipole approximation 

^^i-A = ■ A, (17) 

where Ei is the electric field and Di denotes the electric 
dipole operator for the i-th atom. The electric field con- 
sists of the quantum light field coupling states |0) and |e), 



ijQ*' = i^QQ-'e^'""*, and the classical field coupling states 

|1) and |e), E^^^ = ^;^^^eT^"ot±'fco^ The effective light- 
atom Hamiltonian for the atomic ensemble is obtained by 
summing all the single atom Hamiltonians in Eq. (|17l) . 

iJL-A = Si^^L-A- We introduce the spatially slowly 
varying operators and integrate over space whereby the 
sum over i is absorbed in the operator definitions. Within 
the rotating wave approximation, the final Hamiltonian 
describing the atomic ensemble and its coupling to the 
light fields then reads 



L-A 



A / dVV^O-ee(f,t) - 



■ d^fn{f,t)aei{r,t)-go ^ Bm' 



i)-f H.c. 



mm'7in' 



(18) 



r 



where go = Dq^^-eloquqI c is the coupling constant for 



the quantum light field with Dq — 
is the slowly varying Rabi frequency 



(-) 

e,l 



e^, and fi(r, £) 



(19) 



The first term of the Hamiltonian ([TS]) is the atomic part 
and the second term describes the coupling to the strong 
classical control field. The last term represents the in- 
teraction between the atomic ensemble and the quantum 
light. We see from Eq. (IT5|) that the collective enhance- 
ment leads to an increase of the coupling constant be- 
tween atoms and light by a factor of ^Jn^ ~ V-^A up to 
50 — g^/no, where g is the coupling corresponding to a 
single atom. 

The geometry of the ensemble will lead to coupling 
between different modes mn and m'n' , which is described 
by the coupling matrix 



mn.m'n 



1 

no 



'r,' = — I d'^r±u'^,^{f±)ujn,n'{r±)n{fj_). (20) 



In general, the coupling matrix B„ 



will be z de- 



pendent, but here we have avoided this by choosing a 
basis set fk,mn{r) [Eq. ([2])] in which one can factorize the 
dependence on the transverse f± and the simple exponen- 
tial dependence on the longitudinal position z. Moreover, 
as mentioned earlier, we have chosen the atomic density 
distribution to be constant in the z direction, which leads 
to a z independent coupling matrix Bmn,m'n'- This sig- 
nificantly reduces the numerical complexity in the calcu- 
lations below. In addition, the rotational symmetry im- 
plies that there is no coupling for different values of the 
azimuthal quantum number m, Bmn,m'n' 



Thus, we have an independent set of equations of mo- 
tion for each m value. Note that in the one-dimensional 
limit, where <7± — ?> 00, the atomic density distribution is 
constant riiD = tiq and then the coupling matrix B^} 
transforms to a unit matrix due to the orthogonality of 
the Umn{r±) modes, hence uncoupled modes. 



We assume for simplicity, that the classical driving field 
is spatially uniform with Rabi frequency fl{f,t) = D,{t). 
It was shown in Ref. [l^, that the control field should 
be loosely focused with beam waist at least more than 
twice that of the input light mode, which agrees with our 
assumption of the plane wave control field. Otherwise, 
the position dependent control field will lead to a position 
dependent AC Stark shift, which may lead to reduced 
efficiency of the light storage. 

The above Hamiltonian allows us to derive the 
equations of motion for describing the propagation of 
light through the atomic ensemble as well as the evo- 
lution of the atomic coherences. We use the Heisenberg 
equation and calculate the commutators to obtain 



at dz 



iz,t) 



^mn (^) ^) 



2/co 

ig^c'^B':^] Pmn'{z,t) 

n' 

1 



(21) 



dt 



Prnniz,t) = - -7 + i A P™„(z,t) (22) 



+ -in{t)Smniz,t) 



— Smniz,t) = -ifl*{t)P,n„iz,t). 



(23) 



The procedure employed here applies to a collection 
of near paraxial modes. In reality one should include all 
light modes. The non-paraxial modes would be needed 
if one were to account for the evolution of all sponta- 
neously emitted photons, whereas here we will only be 
trackingthe collectively enhanced emission. It was shown 
in Ref. Q that for a dilute atomic ensemble in the ideal 
gas approximation the effect of non-paraxial modes can 
be accounted for by the inclusion of the polarization de- 
cay due to spontaneous emission at a rate 7, which has 
already been included in the expressions above. We as- 



5 



sume that there is no spin-wave decay that leads to a 
loss of the spin-wave coherence S{z,t). It was shown 
in Ref . [l^l , that such a decay introduces a simple expo- 
nential decay into the solutions for storage and read-out 
and does not make the optimization harder. We also 
omit the quantum noise. In Ref. [23j, it was argued that 
under reasonable experimental conditions, it is not nec- 
essary for calculating any normally ordered product such 
as the efficiency. Alternatively, if required the noise can 
be reintroduced at a later stage. Thus the only two loss 
channels present in our model are the spontaneous de- 
cay 7 of the excited state in other directions than the 
collectively enhanced one and the leakage of the input 
light mode out of the atomic ensemble. The latter one 
results from the fact that not all the light is stored within 
the atomic ensemble, since part of the light will expand 
beyond the atomic ensemble z > L during the storage 
process and hence be lost. For instance, avoiding this 
leakage in the resonant EIT case requires reduction of 
the group velocity Vg ^ i}'^L/{dj) which, however, in- 
creases the spontaneous decay. As a consequence the 
optimal memory performance will be a trade-off between 
these two loss sources. 

Next, we transform the equations of motion into a 
co-moving frame with t' = t — z/c and introduce a di- 
mensionless time t = jt', detuning A = A/7, position 
z = z/L, and Rabi frequency n{t) — Vl{t)/^. In ad- 
dition, we introduce two dimensionless parameters de- 
scribing the properties of the sample: the peak optical 
depth 



do 



7 



and the Fresnel number of the atomic ensemble 



(24) 



(25) 



where Ao is the wavelength of quantum light. The opti- 
cal depth quantifies the absorption of resonant light e^"^ 
and thereby describes the strength of the coupling be- 
tween the light and atomic ensemble. The Fresnel num- 
ber characterizes the geometry of the atomic ensemble 
with regard to the diffraction of light. With this we ob- 
tain the dimensionless coupled equations describing the 
atom-light system 



j~ ^mn (^: ^) 

dz 



,~Pmn {z J i) 

dt 



47rF 

— i\/d.r 

2 



(~z,t) 

(rn) 



(26) 



-^^^/d^^BZlPn^n■{z.t), 



- ( i + i A 1 F„ 



{z,i) + -in{i)S,nn{z,i) 



d 



dt 



]^in*{i)Pran{i,i). 



(27) 
(28) 



Here, a factor of 1/^/7 has been absorbed into amn{z,i) 
and a factor of y/L into the definitions of Pmn{z,i) and 
Smn{z,t). One can see from the above equations of mo- 
tion that the only two physical parameters that describe 
the quantum interface between light and matter are the 
peak optical depth do and the Fresnel number F. Since 
the perpendicular wave number scales as k±_rnn ~ l/ci 
so that the first term in Eq. ([^5]) only depends on the 
Fresnel number F. It will be shown in Section IIIII that 
the parameters of the control field (its temporal shape 
f2(t) and detuning A) do not influence the efficiencies of 
the light retrieval. 

In our model, we neglect the fact that the sponta- 
neously emitted photons can be rescattered and inter- 
fere with the quantum light signal. However, this effect 
is negligible in the paraxial regime, where the solid an- 
gle of emitted light is small. Moreover, the storage we 
are interested in employs only a single or a few photon 
pulses, thus in this case there will be at most a few spon- 
taneously emitted photons. In addition, the probability 
of the spontaneous emission process is small in most ex- 
periments and decreases with increasing optical depth. 

Before proceeding, we can relate the three-dimensional 
theory to the one-dimensional one (2^, where one as- 
sumes that the atomic density distribution is independent 
of the transverse coordinate n(r) = 'n{z) [for our choice 
of the atomic density it corresponds to n{r) = no], which 
leads to an infinite value of the Fresnel number of the 
atomic ensemble F ^ 00. The first term in the propa- 
gation equation pS)) gives either a phase depending on 
the momentum of the incoming quantum light signal or 
it vanishes if the waist of the incoming beam is wide 
enough. In the ID limit, the coupling matrix B^} in 
Eq. (pni) becomes the identity matrix, since the modes 
Umnir'±) are orthogonal. In consequence, the last terms 
in Eqs. ([26]) and ([27]) transform into ^i\/doPmn{z,i) and 
^iy/doamn{z,i), respectively. Since there is no longer 
any coupling between distinct transverse modes, we treat 
each mode independently. In this manner, we end up 
with the one-dimensional equations from Ref. [23| , where 
the only important parameter is the optical depth do. 
The typically employed one-dimensional description thus 
applies to all transverse modes in the limit of _F — )• 00. 



III. OPTIMAL READ-OUT FROM THE 
MEMORY 

In this section, we study the optimal retrieval of light 
stored in a quantum memory based on a A-type atomic 
ensemble. The separate analysis of the read-out process 
is not only interesting for the quantum memory purposes 
but also for gaining knowledge about quantum states of 
atoms generated by other means. Such atomic states can 
be mapped onto light states in order to study their prop- 
erties or quantum correlations [1^. As noted above, the 
equations of motion describing the quantum light-matter 
interface [Eqs. (P5)) - (P5| ] depend only on two physical pa- 
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rameters of the atomic ensemble, F and cIq. In addition, 
the equations of motion also depend on the coupling to 
the control light field r2(t) and its detuning A from the 
transition energy. In the one-dimensional theory, the re- 
trieval efficiency is independent of the detuning and tem- 
poral shape of the control field. In this section, we will 
first show that this holds also in the case of many coupled 
transverse modes. Next, we study the dependence of the 
read-out efficiency on the Fresnel number of the atomic 
ensemble F. 

Initially all atoms are in the ground state and no 
atomic excitations are present, which corresponds to an 
empty quantum memory. The photon storage is based 
on the mapping of the incoming quantum light field onto 
some spin-wave mode in the atomic ensemble. At a later 
time, the stationary excitation is restored onto an outgo- 
ing light mode. The retrieval efficiency is then the ratio 
between the number of retrieved photons and the num- 
ber of stored excitations. If one normalizes the stored 
spin-wave according to J2mnIo ^■S|'5'mn(5, t = 0)^ = 1, 
the read-out efficiency is 



E 



dtK*(i)|2 



(29) 

where we have assumed that the retrieval starts at time 
t = and continues until all excitations are read out from 
the memory at t — co. 

We want to prove that the read-out efficiency r],- de- 
pends neither on the detuning A nor on the temporal 
shape of the control field To this end we use the 

conservation law 



dz dt dt 

(30) 

which states that the number of excitations in the system 
can change only due to the spontaneous emission, since 
this is the only decoherence channel for the read-out pro- 
cess. This can be derived directly from Eqs. 
Here, we have for simplicity introduced a compact vec- 
tor notation '*™"('^' ^) ~^ d^™\z,i), which replaces 
the mode notation with indices and sums. Next, we in- 
tegrate Eq. ([30|) over space and time ^ dz dPj . We 
use the initial conditions that there is a normalized spin- 
wave excitation in the ensemble S{z,t = 0) = S'o(5), 
/p dz\SQ{z)\^ = 1 and no coherence between the ground 

and excited states P(z, 0) = 0, as well as the boundary 
condition that there is no incoming quantum light field 
during the read-out procedure: din{z = 0, f) = 0. We also 
assume that we deal with a complete retrieval, where no 
excitations are left in the system after the read-out is 
finished: S{z, oo) — P(z,oo) = 0, which can be ensured 
by sending in sufficiently strong control light. Then we 
find that the retrieval efficiency also can be expressed as 



unity minus the losses during the procedure 

/• 1 (*OQ 

r/r = 1 - / dz I di\P{S,i)\^ 



(31) 



= 1 - 



dz l{z), 



where l{z) is a position-dependent loss per unit length: 



di\P{S,i)(' 



(32) 



£-1 



diP\u\i)P{u,i) 



and C ^ denotes double inverse Laplace transform 
{u — >■ z,u' ^ z') in the position coordinate, f{u) — 
C[f{z)] — dze~^^f{z), evaluated at position z in both 
cases. In the following part of the paper we will use 
the purely imaginary u and transform to u — > iu, where 
from now on u is purely real. Next, we wish to express 
the loss in terms of the Laplace transformed spin-wave 
So{u) = S{u, t = 0) by a relation of the form 



1{S) = C-^ Sl{u')A{u',u)S{u) 



(33) 



In view of the boundary conditions, this can be achieved 
if we can find a matrix A(w, u') that fulfills 



d r 
di 



PHu', i)A{u', u)P{u, i) + S^u', i)A{u', u)S{u, i) 

= P\u',i)P{u,i). (34) 

The compact notation of the matrix A{u',u) is short- 
hand for a four-dimensional matrix with components 
^'Z'iu',u). InAppendixEl we prove that it is indeed al- 
ways possible to find a matrix A{u' , u) fulfilling Eq. (|34p . 
This allows us to rewrite the position dependent loss in 
terms of the spin-wave So{u) and the matrix A(u,u'), 
which only depends on the optical depth do and the Fres- 
nel number F while independent of the detuning A and 
the temporal shape of the driving field f2(t). The loss is 
proportional to the ratio between the decay of the spin- 
wave into the unwanted directions and the constructively 
interfering decay into the desired output light pulse with 
a direction set by that of the control beam. Note that in 
the case of nonzero spin dephasing 7^, which introduces 
exponential decay of the spin- wave and the corresponding 
output light, the loss becomes dependent on the control 
parameters and the shape of the output pulse [23| . 

After finding the matrix A{u, u') as the solution of the 
Sylvester equation (|Aip . we can rewrite the expression 
for the read-out efficiency (j3ip as a matrix product cor- 
responding to an eigenvalue problem 



'7r 



A{u' ,u)]Sq{u)^ (35) 



du 



du'Sliu') 



AT:^i{u — u') 



X [I- A{u',u)]Sq{u) 
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FIG. 2. (Color online) Maximum retrieval efficiency as a func- 
tion of the Fresnel number of the atomic ensemble F, for 
different transverse azimuthal quantum numbers m, and for 
fixed peak optical depth do = 40 (top) and do = 100 (bot- 
tom). The dashed black lines denote the ID limit of _F — > oo. 



In order to calculate the efficiency numerically, we dis- 
cretize the Laplace transformed position z — ^ u with a 
finite cut-off Umax- The value of this cut-off and how fine 
the grid resolution needs to be depend both on the phys- 
ical parameters and F . A finite number of Bessel 
modes depending on the Fresnel number F has been 
used in order to set the coupling matrix B in Eq. (PD)) 
in the Sylvester equation (|A1|) . This equation is solved 
numerically to find the matrix A(u', u). The integral over 
the position z can be calculated analytically. Next, we 
construct the matrix resulting from the inverse Laplace 
transform [e^^""''') - l]/[AT:'^i{u - u')]. The two inte- 
grals over u and u' are replaced by sums and we re-index 
the four-dimensional matrix A so as to turn each pair 
of indices into one index A„„'(u,u') — > An^n', where 
n = {n,u) and h' = {n' ,u'). Thus, the sums of the 
four-dimensional matrices can be implemented as mul- 
tiplications of two-dimensional matrices. The retrieval 
efficiency is thus an eigenvalue of the large kernel ma- 
trix - [I - A(w',m)] [e'("-"') - 1\/[At:H{u - u')] with 
the corresponding eigenstate Sq{u). The optimized spin- 
wave to retrieve from is the eigenvector corresponding 



to the highest eigenvalue. In turn, we can determine 
the corresponding optimal outcoming light mode. An 
alternative method of calculating the retrieval efficiency 
(described in detail in the next section) involves the adia- 
batic elimination of the excited state followed by calcula- 
tions with the Laplace transform in the time variable t. It 
is worth mentioning that this optimal spin-wave to read- 
out from is not necessarily the best one for the combined 
process of write-in and read-out. Thus the full memory 
procedure has to optimized, which is done in detail in 
SecHVl 

In the case of the one-dimensional theory, the read- 
out efficiency depends only on one physical parameter: 
the optical depth do as ?7r ^ 1 ^ l/c^o in the dp — oo 
limit. Here, we include an additional degree of freedom 
that describes the geometry of the atomic ensemble with 
regards to the diffraction of light, namely the Fresnel 
number of the atomic ensemble F. 

We have optimized the read-out process from the mem- 
ory numerically for two fixed values of the peak optical 
depth do = 40 and do — 100 and we have analyzed the 
influence of the Fresnel number F . Here, the memory 
operates on resonance, A = 0. The maximum retrieval 
efficiency for different branches of the azimuthal quan- 
tum number m is shown in Fig. [5] for do — 40 (top) and 
do = 100 (bottom). The results are plotted together 
with the one-dimensional case of infinitively large Fres- 
nel number, F ^ oo (dashed black lines). The maximal 
retrieval efficiency of the spatial memory approaches this 
limit already for small F ^ 1 for several values of m. 
The figure clearly shows the spatial multimode character 
of the memory, since high efficiencies for different values 
of the azimuthal quantum number m exist. Furthermore, 
in the plot we show only the optimal mode for each m. 
In general, there are several orthogonal modes with high 
efficiencies within each m-subspace. 

In Fig. [3l we show the density of the optimal spin- 
wave excitation to read-out from |5opt(/5, 5)p and the 



corresponding optimal output light intensity 



where p ~ pja ^. Here, we show the results for the Fres- 
nel number F — 0.5 and for low (do = 0.1) and high 
(do — 100) values of the peak optical depth. For sim- 
plicity, we assumed a constant Rabi frequency of the 
classical light field fi. Note that the radial position is 
scaled with respect to the width of the atomic ensemble 
p = p/<y± and time is dimensionless and scaled according 
to t — t|r2p/7. In the case of a low absorbing medium, 
do = 0.1 [Fig. [31 (left)], the spin-wave is spread over the 
entire ensemble length and its density is symmetric with 
respect to the middle of the sample in the axial posi- 
tion z. For high optical depth do = 100 [Fig. [3l (right)], 
most of the excitation is localized closer the end of the 
sample {z = 1) from which the light is read out. In the 
latter case, the excitation does not need to propagate 
through the entire length of the sample, which in conse- 
quence leads to lower losses due to spontaneous emission. 
Thus, if the optical depth is high enough it is better to 
effectively use only a part of the ensemble and thereby 
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FIG. 3. (Color online) (top) The intensity of the optimal 
output light modes for retrieval at the end of the atomic en- 
semble, z = I for (left) low do = 0.1 and (right) high do ~ 100 
peak optical depths, (bottom) The corresponding density of 
the optimal spin-wave excitations stored in the atomic mem- 
ory. In both cases, the Fresnel number of the ensemble is 
F — 0.5 and retrieval is resonant A = 0. The radial position 
is scaled with the atomic ensemble width as p = p/^± and 
time according to i — t\Q,\^ /■y. 
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FIG. 4. (Color online) Schematic illustration of the two op- 
eration modes of a quantum memory: (top) forward, where 
the control fields il{t) during the storage and read-out are in 
the same direction, and (bottom) backward read-out with the 
opposite directions of the control fields. 



ical parameters: the peak optical depth do and the Fres- 
nel number F. In order to understand the operation of 
a full quantum memory, one also needs to store the light 
beforehand and we now turn to the full process of stor- 
age followed by retrieval. A quantum memory is often 
discussed in two different modes of operation (see FigH]): 
forward operation, where the stored light is retrieved in 
the same direction as it was written in and backward op- 
eration, where the read-out is performed in the opposite 
direction [ll.l23t. 



Read-out in the forward direction 



decrease the propagation length. The distribution of the 
optimal spin-wave excitation in the transverse direction 
reflects the atomic distribution, but also depends on the 
Fresnel number of the ensemble, which affects the diffrac- 
tion of the input and output light pulses (see detailed 
discussion in Section HVB I) . 

The transverse shape of the corresponding optimal out- 
put light pulse also reflects the Gaussian distribution of 
the atoms. However, its duration and temporal shape 
strongly depend on the peak optical depth do , which was 
studied in detail in the ID theory 23|. For the resonant 
case under consideration, A = 0, the output pulse length 
is proportional to the optical depth do: since it is propor- 
tional to the inverse of the group velocity tout ~ ~ 
(io7/|l^p. Note that for the off-resonant Raman transi- 
tion (A ^ 7^0), the pulse duration dependence on the 
optical depth is different, tout ^ ^^/('^o7l^P)- Since the 
retrieval efficiency is detuning independent, the choice of 
the experimental realization may be affected by the avail- 
able laser power and the length of the pulses one wants 
to retrieve. Unlike the efficiency, the shape of the light 
pulse does depend on the detuning A. 



IV. ADIABATIC STORAGE AND RETRIEVAL 

In section IIIIl we solved the full three-dimensional 
problem for the read-out of a spin-wave from an atomic 
ensemble and studied its dependence on the crucial phys- 



In this subsection, we derive analytical expressions for 
the efficiency of storage followed by retrieval in the for- 
ward direction together with the input-output relations 
between the incoming and outcoming light modes and 
the spin- wave in the atomic ensemble. 

For simplicity, we restrict ourself to the limit where we 
can adiabatically eliminate the excited state by setting 
■^P{z,i) « 0. The condition for this is ipuiserfo7 ^ li 
where tpuiso is the duration of the preferably smooth 
quantum light pulse. The adiabatic elimination leads to 
the following equations of motion (now written in com- 
pact vector/matrix notation) 



P{~z,i) 



-T:a{z,t) 
dz 



i + zA 



M{S,i), (36) 



AttF 
/d^n{i) 



iA 



d{z, t) 



(37) 



-~-S{z,t) 



i + zA i + zA 



5a(z,ij(38) 



Here, the imaginary part of the first term in Eq. (j37p 
is the phase shift due to the index of refraction of the 
medium with rzrcfr = - i k'^^a / {4n F) + do AM'^/{1+ 4 A'^). 
Note that fc^ although written as a vector is a diago- 
nal matrix with /c^ values on the diagonal. The AC 
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Stark shift of the atoms is accounted for by the imag- 
inary part of the first term in Eg. (1551). A\n(t)\^/(1 + 
4A^) and is zero in the resonant case. The real parts 
of these two terms describe the damping of the hght 
and spin-wave excitations with the respective absorp- 
tion coefficient ^doM'^/(l + 4A^) and effective decay rate 

-I- 4A^) due to the spontaneous emission at 
a rate 7. The last terms in Eqs. (|57)) and ([551) represent 
the coherent interaction between light and the atomic 
ensemble. 

It is worth noting that the control field deter- 
mines the speed of the memory and can be completely 
eliminated from the equations of motion by the rescal- 
ing: P{z,i) P{z,t)IVt{i) and a{z,i) a{z,i)/n(t) 

and changing variables i ^ v — h{0,i) ~ di'\tl{i')\^ , 
where h{0, t) is the integrated intensity of the control 
classical field. Thus, the dynamics of the system is inde- 
pendent of the temporal shape of the driving light field 
described by which merely reshapes the shape of 

the incoming light mode which is stored and accordingly 
the shape of the output mode. However, here we do not 
perform the rescaling in order to avoid introducing addi- 
tional notation. 

From the equations of motion ([57]) and (1551) . one can 
derive expressions for the spin-wave excitation stored in 
the sample and for the output light for given incoming 
quantum light pulse. We have solved these equations in 
two different ways. In the first case, we use the Laplace 
transform of the position z and we obtain equations in 
momentum space u and time t. This transformation was 
also used in the previous section for proving that the 
read-out efficiency only depends on the optical depth dp 
and the Fresnel number F. The second approach is based 
on Laplace transform in time and leads to solutions in po- 
sition space and frequency. The numerical calculations 
were performed with both approaches yielding the same 
results. Here, we only discuss in detail the analytical 
expressions for the forward read-out with Laplace trans- 
form in space. The alternative approach as well as the 
backward mode of operation are described in detail in 
Appendicies |B] and O respectively. 

Solving the coupled equations of motion ([M]). ([Tf]). 
and ([55]) . the connection between outcoming light, spin- 
wave, and incoming light can be written in the form of 
input-output beam splitter relations. In addition to the 
initial and boundary conditions introduced for the read- 
out in Sec. IIIIl we use here the appropriate conditions 
for the storage procedure. There is no spin-wave excita- 
tions in the ensemble at the start of the writing process 
S{z, t = —00) = and no coherence between the ground 
and excited state P{z,t = —00) = 0. The storage of the 
incoming light field d{z — 0,t) — ain(i) takes place in the 
time interval t € (— oo,0] and is normalized according 
to J_^dt\din{t)\'^ — 1. We perform Laplace transform 
in space z, which leads to elimination of the differential 
equation ([57)) for the quantum light d{z,t). Thus, the 
above mentioned relations in momentum and time space 



read 



flout (?) = ^ 



So{u) 



dt 



du M[h{i), i, u]e™ Soiu), (39) 



^[n*{-i),-i,u]din{i). (40) 



Here, M[n{i),i,u] = Q(u)eN(")* with 



N{u) 
T{u) 



i+zA 



(41) 
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do|f^(?)|' 



zA 



fc2 ^2 



AttF 



(42) 



and the coupling matrix B is defined in Eq. ([^0)) . 

We start by considering the process of writing infor- 
mation into the atomic ensemble. The storage efficiency 
is the ratio between the number of stored excitations and 
the number of incoming photons. If the incoming light is 
normalized to one, the storage efficiency is given by 



Vs 



Z\\2 



dz\Soiz)\ 



(43) 



du 



du' 



A'K'^i(u' — u) 



Sliu)So{u'). 



We substitute the corresponding input-output relation 
from Eq. ([40]) and obtain 



Vs 



du / du' / di dt' 



(44) 



[n*{-t),~i,u] 



ATr^i{u' — u) 

xMF[n*{-i'),~i', u']din{i'). 

In order to numerically calculate the efficiencies, 
we discretize time t and momentum u and introduce 
finite cut-offs in both of them. The requirements for 
the cut-off values depend on the chosen peak optical 
depth do and the Fresnel number F. Next, we construct 
the M matrix and the term resulting from the inverse 
Laplace transform. The sizes of these depend on the 
number of discrete t and u points we take. Now, the 
integrals can be replaced by sums and, in consequence, 
can be implemented as matrix multiplications. In 
the end, we have a large kernel matrix Cs[t^t'] ^ 
M*[ii*(-f),-?,u]M'^[ii*(-?),-t',w'][e*("'-") 
l]/[47r^i(u' — u)], which after reindexation can be 
diagonalized in order to optimize the storage. The eigen- 
values of this matrix give the storage efficiencies with 
corresponding incoming quantum light modes ain(i)- 
The highest eigenvalue is the maximal storage efficiency 
and the corresponding eigenvector is the optimized 
incoming light mode for storage. 
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Now, we turn to the full combined procedure of storage 
of an incoming light pulse followed by forward retrieval 
of the spin- wave. The relation between the incoming and 
outcoming light at z = 1 is 



flout (?) = 



1 



du I di' M[fi(i),i,M] (45) 



2TTi 

xWf[n*{~i'),~i',u] din{i') 



The efficiency of this forward quantum memory is thus 
given by 



30 /-OO />0 />0 

du / du' / di di' 

OO J — OO J — OO J —OO 



1 . 

xalii) M*[n*{~i),-i, u] Mi[n{i^), i, (46) 
xM^[n{i,),i,,u'] M^[n*{-i'), -i', u']din{i'), 

where M„[f2(t), i, u] = M[f^(t), f, u]e'". Analogously to 
calculating the efficiencies for storage only, as discussed 
above, we discretize time and momentum in order to 
allow for numerical diagonalization of the large kernel 
matrix Cs+h[i,i'] ~ M*[r2*(-t), -f, w] Mlp {i^),ir,,u] 
M„[ii(i),4,u'] U^[n*{-i'),^i',u']/{4:n^). This gives 
the efficiencies for the forward memory with correspond- 
ing incoming light modes ain(i)- 

An alternative approach is to use Laplace transform in 
time and derive the beam splitter relations in the posi- 
tion and frequency domain. This is described in detail 
in Appendix [Bj The analytical expression for the back- 
ward read-out for both solution methods are presented 
in Appendix ICl 



B. Results for backward and forward memories 
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FIG. 5. (Color online) The maximal efficiencies for storage 
followed by read-out in the forward (top) and backward (bot- 
tom) direction as a function of the Fresnel number F for dif- 
ferent values of the azimuthal quantum number m and for a 
peak optical depth do = 100. The dashed black lines denote 
the corresponding ID limits with F — >■ oo. 



In the previous subsection as well as in Appendices |B] 
and \C\ we derived the input-output beam splitter rela- 
tions between the spin-wave and the incoming and out- 
coming light modes as well as expressions for the efficien- 
cies for both operation modes (forward and backward). 
In this subsection, we present the results of numerical 
calculations of the maximal efficiencies and the optimized 
light modes and spin-waves. In the one-dimensional the- 
ory, the optimal storage is the time reversal of the optimal 
retrieval. Since the optimal spin- wave excitation is real in 
that case, the combined process of storage and retrieval is 
independent of detuning A [2^ . In the three-dimensional 
theory, however, the spin- wave is not real and one cannot 
directly use the time reversal argument, which will result 
in dependence on the detuning A. Here, for simplicity 
we only consider the resonant case with A = 0. 

In Fig.[Sl we show the maximal efficiencies for resonant 
(A — 0) memories operating in the forward (top) and 
backward (bottom) direction as a function of the Fresnel 
number F, where the peak optical depth has been fixed 
to do = 100. The highest memory efficiencies for differ- 
ent azimuthal quantum numbers m = 0, 1,2, and 3 are 
shown. In addition, for each value of m, there are several 



orthogonal modes with high efficiencies which are not 
shown in these plots. For both directions of the read- 
out, the simulations show the multimode character of 
the memory, i.e. there are many orthogonal modes even 
within each m-subspace to store into and retrieve from 
with high efficiency. Note that the negative values of 
the azimuthal quantum number m have the same effi- 
ciencies as |m|. The efficiencies decrease for increasing m 
due to the spatial light distribution in high m modes. 
Light beams with |m| > vanish towards the center of 
the sample (p — >■ 0) as p'™'. For large m, it is there- 
fore harder to focus the input light into the dense center 
of the atomic ensemble, which results in lower effective 
optical depth. Nevertheless, for large Fresnel numbers 
F > 1, the efficiencies for successive m decrease only by 
a few percent in comparison to m — 0. For the Fres- 
nel number F = 6, the efficiency difference between the 
modes with subsequent m values is smaller in the case 
of forward read-out AT^for ~ 0.1% in comparison to the 
backward one with Ar^back ~ 0.6%. We compare the re- 
sults of the three-dimensional theory with the efficiencies 
obtained for the one-dimensional memory (dashed black 
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FIG. 6. (Color online) The maximal efficiencies for storage 
followed by backward (full circles) and forward (empty cir- 
cles) read-out as functions of the Fresnel number F for optical 
depth do = 40 (bottom, dotted green lines), do = 100 (middle, 
dashed-dotted red lines), and do = 200 (top, solid blue lines). 
The dashed lines for large values of the Fresnel number F de- 
note the corresponding values from the one-dimensional the- 
ory. 



lines). The efficiencies approach the ID limit of F — )■ cxd 
already for small Fresnel numbers _F w 1. Moreover, the 
one-dimensional limit is reached faster for the forward 
memory. 

In Fig. [6l we compare the highest achievable efficien- 
cies for memories operating in backward (full circles) and 
forward (open circles) directions. The efficiencies are 
plotted as functions of the Fresnel number F for three 
values of the peak optical depth do — 40, 100, and 200. 
The dashed lines for large Fresnel numbers F denote 
the corresponding values from the ID limit. The one- 
dimensional theory predicts that the backward memory 
always yields higher efficiencies [2^. Here, we find that 
forward read-out is advantageous for small Fresnel num- 
bers. The intersection between the backward and forward 
efficiencies appears at lower Fresnel numbers for increas- 
ing values of the peak optical depth d^. The advantage of 
forward read-out for small Fresnel numbers F can be un- 
derstood by considering the diffraction of the light beam 
(in analogy to Gaussian beam diffraction) . In the case of 
backward retrieval, the z dependent part of the phase of 
the excitation is reversed by means of the classical control 
field with propagation direction opposite to the incoming 
quantum field. However, since we are not able to phys- 
ically complex conjugate the stored spin-wave, the evo- 
lution of the transverse phase profile will be unchanged 
and will lead to reduced constructive interference of the 
restored spin-wave: the output light will accumulate an 
unwanted spatially varying phase, which in the case of 
a Gaussian beam can be written as ^ exp [ik^r^^/ R{z)]. 



Here, R{z) is the radius of curvature which is propor- 
tional to the Rayleigh range Zq as R{z) ^ Zq. In the 
regime of small Fresnel numbers, the far field divergence 
angle 8 ^ ^JF\/L is small and the resulting Rayleigh 
range zo ^ ttLF is short. Thus for small Fresnel numbers 
where the diffraction is large, the unwanted phase grows 
since it scales as ~ 1/F^. The desired retrieval process 
is therefore not phase matched and cannot conserve mo- 
mentum. This reduces the efficiency of the storage fol- 
lowed by backward retrieval. In contrast, for the forward 
retrieval we do not have this problem with the transverse 
phase profile and thus forward read-out is less sensitive 
to small Fresnel numbers F. This discussed reduction 
of the backward retrieval efficiency is analogous to the 
reduction due to accumulated phase for non degenerate 
ground states found in Ref. [2^. In the case of large 
Fresnel numbers which is well described by the ID the- 
ory, the profile of the excitation is flat in the transverse 
direction and no unwanted phase is accumulated. Thus 
for large Fresnel numbers, backward retrieval is favor- 
able. For high peak optical depth do it is feasible to use 
only part of the atomic ensemble, whereby the effective 
Fresnel number is increased. Therefore the intersection 
between forward and backward operating memories for 
higher peak optical depths do occurs at lower values of 
the Fresnel number F. 

It is interesting to note that in the case of the DLCZ 
protocol, where the spin- wave is generated by a paramet- 
ric gain interaction and next released from the ensemble 
with the beam splitter interaction, the above described 
mode mismatch (via the unwanted spatially dependent 
phase accumulation) will occur for the read-out in the 
forward direction ^30.] . In that case, the backward mem- 
ory will work better than the forward one for small values 
of the Fresnel number F. 

It was shown in the ID theory that in the limit of 
high optical depths do ^ oo the inefficiency scales as 1 — 
771D ~ 1/do- We have fitted the minimal inefficiencies 
in the limit of high values of the Fresnel number F < 1 
with the power dependence on the Fresnel number 1 — 
V3D ~ (l-r7iD)(l + l/i^') ~ (l + l/i^O/4- For backward 
operating memories, we obtain 1 = 3/4 and for forward 
retrieval I = 0.9. 

Next, we fix the Fresnel number and consider the de- 
pendence on the peak optical depth do- In Fig- 13 we show 
the maximal efficiencies for the backward (black) and for- 
ward (light blue/gray) memories for two values of the 
Fresnel number = 0.2 (open circles) and F = 2 (stars) 
together with the corresponding results in the ID limit 
(dashed lines). For the smaller Fresnel number F = 0.2 
one can see that for low optical depths do the efficien- 
cies are considerably lower than the ones corresponding 
to the ID limit, which in the case of backward read-out 
is approached for high do « 300. The efficiencies for the 
larger Fresnel number F = 2 approach the ID limit al- 
ready for small optical depths especially for the forward 
operating memory. The backward and forward efficien- 
cies for F = 0.2 intersect at the optical depth do = 10 
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FIG. 7. (Color online) The maximal efficiencies for the back- 
ward (black) and forward (light blue/gray) memory as a func- 
tion of the peak optical depth do for F = 0.2 (empty circles) 
and F — 2 (stars). The dashed lines correspond to the results 
of the ID theory with F — ^ oo. 



whereas for larger do, the backward memory is advanta- 
geous. 

The three-dimensional theory presented here allows for 
calculation of not only the efficiencies but also the opti- 
mal light modes to store into and retrieve from the mem- 
ory. The optimal input and output light modes can be 
shown to be time reversed complex conjugates of each 
other aout{t) ~ a*^{—t), which we prove in Appendix IdI 
This property can be used for an iterative procedure to 
obtain the maximum efficiency of the quantum memory. 
First, one stores an arbitrary input light mode and re- 
trieve it. In the next steps one uses the time reversed 
and complex conjugated output light mode as an input 
mode for the next run of the experiment. Such successive 
time-reversal iterations with complex conjugation lead to 
optimal light storage. This is analogous to the proposed 
procedure in the one-dimensional case [23|, which has 
been demonstrated experimentally in a warm Rubidium 
vapor [11]. 

The light modes will, in general, be correlated in the 
spatial and temporal degrees of freedom and the spin- 
waves in the transverse and longitudinal spatial modes. 
This means, for instance, that the transverse spatial pro- 
file of the optimal incoming light mode would change in 
time. Such correlated light modes are disadvantageous 
from an experimental point of view, where it is simpler 
to have a time independent transverse mode profile. 

In order to investigate the degree of correlation for 
the optimized light modes we have calculated their puri- 
ties in terms of the Schmidt decomposition. Any mode 
can be factorized as 0;°''* ("tx , ?) = J2n hn{r±)an{i), 
where hn{r_\_) and a„(t) depend exclusively on the trans- 
verse component and time, respectively, and Pyi — 1. 
The largest P„ is a measure of the purity of the mode. 
In the ideal pure case, there is only a single compo- 
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FIG. 8. (Color online) (top) Purity of the optimal input light 
mode for backward (solid lines) and forward (dashed lines) 
memory as a function of the Fresnel number F for the peak 
optical depth do = 100 and different values of the azimuthal 
quantum number m = Q (circles), m = 1 (squares), m = 
2 (triangles), and m = 3 (stars), (bottom) The purity of 
the optimal input light mode for backward (solid lines) and 
forward (dashed lines) memory as a function of the Fresnel 
number F for m = and different values of the peak optical 
depth do = 40 (circles), do = 100 (squares), and do = 200 
(triangles) . 



nent in this expansion which consequently has P = 
1. In order to determine the purity of the optimized 
modes, we construct the reduced density matrix g{i' , i) = 

C'''^(^^i-,^>ir(^^i-,*0/Tr[ar^(r='^,t>r(^^^,i')] the 
eigenvalues of which are exactly the set Anal- 
ogously, one can calculate the purity for the stored 
spin-wave excitation which can be decomposed as 

^opt(rL, Z) = Y.n \fPli fn{r±)Sn{S). 

In Fig. |5] (top) , the purity of the optimal input light 
mode for various azimuthal quantum numbers m is shown 
as a function of the Fresnel number F for a peak optical 
depth do — 100 for both forward (dashed lines) and back- 
ward (open lines) operating quantum memories. The pu- 
rities of the optimal output and input light modes are 
equal since the output light is the time reversed complex 
conjugate of the input mode (see Appendix iPj) . First 
of all, the purities are in general very close to one. For 
large values of the Fresnel number F, the purity is higher 
in the case of the forward memory even though the spin- 
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FIG. 9. (Color online) (top) The spatial profile of the intensity 
of the optimal incoming light pulse for F = 0.2 (black lines) 
and F = 2 (light blue/gray lines) for backward (dashed lines) 
and forward (solid lines) memory (do = 100). (bottom) The 
corresponding spatial profile of the phases as in (top). 



wave is not factorizable. It is interesting that in all cases, 
there exists a minimum of the purity which depends on 
m. If one decreases the Fresnel number, the purity drops 
at most by a few percent and is bounded from below at 
^ 96%. Next, we investigate how the optical depth in- 
fluences the purity, see Fig. |8] (bottom). For large F, 
the purities have the same value for all three values of 
the peak optical depths do = 40, 100, and 200. The sit- 
uation changes for small where the purity decreases 
more for large values of do. Again, in this regime, the 
purity for the forward read-out decreases to lower val- 
ues than the backward one. All in all, the purities of 
the optimal light pulses are close to one and their depen- 
dence on the physical parameters do and F as well as the 
choice of the azimuthal quantum number m is very weak, 
which is advantageous for experimental realizations. One 
can, however, estimate how the fact that one uses only 
the dominant Schmidt component of the optimal input 
light pulse and measures only the Schmidt component of 
the corresponding output light will affect the efficiencies. 
The efficiency of the full optimized light rj will at most 
be reduced to Typuro ~ '?(! — 4e), where £ = 1 — P is the 
lack of purity (see Appendix |E] for derivation). 

Given the fact that the optimal light modes have very 
high purities, it is reasonable to extract the dominant 



Schmidt component hn{r±_). In Fig. |9] (top), we plot 
the spatial intensity profile of the optimized incoming 
light pulses for small F ~ 0.2 (black lines) and large 
F — 2 (light blue/gray lines) Fresnel numbers, for back- 
ward (dashed lines) and forward (solid lines) memory. 
The transverse spatial shape of the optimal light pulses 
with the azimuthal quantum number m = is well ap- 
proximated by a Gaussian beam. For atomic ensembles 
with large Fresnel number, the optimized light is more lo- 
calized in comparison to the sample size than for small F. 
In the latter case, too narrow pulses would entail a strong 
diffraction leading to excitation leaking out of the sam- 
ple and weaker interaction with atoms. Thus, larger val- 
ues of the Fresnel number allow for more localized light 
pulses, which in consequence will lead to a larger number 
of available modes with high efficiencies. The full multi- 
mode capacity of spatial memories based on the A-type 
atomic ensembles will be studied in detail elsewhere !31| . 
The spatial profiles of the intensity for backward and 
forward memory are only slightly different. The width of 
the pulse is smaller for the backward read-out for small F 
and larger for large F due to diffraction. 

We also show the spatial phase profiles of the optimal 
input light, see Fig. [9] (bottom). The lines correspond to 
the discussed intensities shown in Fig. [3] (top). One can 
see, that for both small and large F, the phase profile for 
backward memory is flatter than the one for the forward 
read-out. It results from the fact, that the focal plane 
with a fiat phase profile of the optimal stationary exci- 
tation for the forward read-out lies in the middle of the 
sample at z = 0.5 while for the backward case it will be 
placed much closer to the read-out end of the ensemble. 
Thus the propagation length is longer for the forward 
retrieval and the accumulated phase is larger. 

The information gained from the spatial profiles of the 
intensity and phase allows us to calculate the beam waist 
of the optimal input light mode and the distance from the 
beginning of the ensemble to the focal plane for the light 
propagating without the atoms. We fit the dominant 
Schmidt component of the optimal mode with a Gaussian 
beam /ipurc(p) ^ exp[—p^/w^{z))]exp[ikp'^/{2R{z))]. 
Here, w{z) — wo\/l + z'^/zq is the dimensionless spot 
size and w{z) — w{z)/a±, wq — wq/(t± is the dimen- 
sionless beam waist, p — p/(7±, and k — kL. The phase 
profile of the input mode is fit to find the dimensionless 
radius of curvature R{z) = z + z'q/z. From these we 
calculate the beam waist 



Wq 



1 



i?2(5f) w^{zf)' 



(47) 



and the distance of the focal plane from the beginning of 
the sample 



Zi ■ 



iJ:^w'^FR{zi) 



(48) 



In order to check the quality of the fitting procedure, 
we calculated the overlap between the pure part of the 
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FIG. 10. (Color online) (top, left) The dimensionless waist of 
the optimal input light multiplied by the square root of the 
Fresnel number F, woVF — ■wo/{\oL), as a function of the 
Fresnel number for backward (solid lines with markers) and 
forward memory (dashed lines). Here, we present the fits for 
the peak optical values of do ~ 40 (blue circles) , do = 100 (red 
squares), and do — 200 (green stars), (top, right) The dimen- 
sionless waist of the optimal input light scaled as in (top, left) 
as a function of the peak optical depth do for forward read- 
out with F — 2 (dashed blue line) and for backward read-out 
with F = 0.2 (solid green line) and F = 2 (dashed-dotted red 
line), (bottom, left) The dimensionless distance between the 
beginning of the atomic ensemble and the focal plane in the 
free space zt = Zi/L of the optimal input light mode as a func- 
tion of the Fresnel number F. (bottom, right) The distance 
Z{ as a function of the peak optical depth do. 



optimal light mode and the fitted mode. These overlaps 
are equal or larger that 0.996 and 0.994 for forward and 
backward operating memory, respectively. 

In Fig. [TU] (top, left), we plot the dimensionless waist 
of the pure part of the input light beam multiplied by the 
square root of the Fresnel number F, iboy/F — ujo/y/ \oL, 
as a function of the Fresnel number. For the forward 
operating memory (dashed lines), we can see a univer- 
sal behavior that the waist is independent of the Fres- 
nel number, it is a constant equal to wq = 0.3\/XoL w 
l/{n^/XoL), given by the wavelength Aq and the length 
of the sample L. For the optimal light beam, it is desir- 
able to both have equal waist along the entire atomic en- 
semble to provide high collective interaction with atoms 
and to be as focused as possible. The trade-off between 
these effects is achieved for an optimal waist which has 
a Rayleigh range zq ~ L independent of F. The sit- 
uation is different for the backward read-out (solid lines 
with markers) , where the spin- wave excitation is not sym- 
metric in the z direction. Here, the optimal beam waist 
grows for increasing Fresnel number. The dependence of 
the beam waist on the peak optical depth is shown in 
Fig. [TU] (top, right). One can see that it is again constant 
for the forward read-out whereas for the backward one it 
decreases slowly with growing optical depth do- 
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FIG. 11. (Color online) (top) The spatial and temporal shape 
of the optimal outcoming light pulse, (middle) The opti- 
mal spin-wave excitation in the atomic ensemble, (bottom) 
The corresponding incoming light mode, (left) The backward 
memory with read-in and read-out at 5 = 0. (right) The for- 
ward memory. All results are for F = 2 and do — 100. The 
radial position is scaled with the atomic ensemble width as 
p = p/cTx and time according to i — t\Q,\^ /j. 



The corresponding distance Zf between the focal plane 
and the beginning of the ensemble is shown in Fig. [10] 
(bottom, left) first as a function of the Fresnel number F. 
For the forward read-out, the focal plane is always in the 
middle of the sample {z = 0.5). In the backward operat- 
ing memory, the focal plane is closer to the beginning of 
the ensemble and decreases only slightly for small Fres- 
nel numbers F < 1. For large values of F it reaches a 
constant value depending on the peak optical depth. The 
influence of the optical depth on the distance Z{ is shown 
in Fig. [To] (bottom, right). Again, for the forward oper- 
ating memory, it is constant and equals to Z{ = 0.5. In 
the case of backward retrieval, high optical depths allow 
for focusing closer to the beginning of the ensemble. The 
distance increases for decreasing peak optical depth do 
as the spin- wave becomes more and more symmetric due 
to the increasing transparency of the medium. 

Knowing the spatial profile of the intensity and phase 
of the optimized incoming light mode, it is also inter- 
esting to study the spatial profiles of the corresponding 
spin-waves as well as the spatio-temporal shapes of the 
outgoing light modes. In Fig. [TT] we show the full set of 



15 



light modes and spin-wave excitations for the backward 
(left column) and forward operating memory (right col- 
umn) with Fresnel number F = 2 and with peak optical 
depth do = 100. We assumed for simplicity a constant 
Rabi frequency and the memory is resonant A = 0. 
In the case of the backward memory, it is preferable 
(especially for high rfo [23]) to store most of the exci- 
tation closer to the beginning of the atomic ensemble 
(^z — 0) where the read-in and read-out are carried out, 
see Fig. [TT] (middle, left). Thus, the light does not need to 
propagate through the entire ensemble, which leads to re- 
duced decoherence due to spontaneous emission as well as 
there being less irreversible phase acquired. Such a spin- 
wave is obtained with an incoming light mode, where the 
time dependence is a "cut/half Gaussian with the in- 
tensity maximum at the end of the write-in process, see 
Fig. [TT] (bottom, left). This is in a gree ment with the re- 
sults of the one-dimensional theory [23] and has been also 
demonstrated experimentally [llj | . For the forward mem- 
ory. Fig. [TT] (middle, right), it is advantageous to store 
a symmetric and smooth spin-wave, which is achieved 
with a light pulse with a Gaussian type time dependence 
(with a very small final cut-off), see Fig. [TT] (bottom, 
right). As mentioned earlier, the outgoing light modes 
are the time reversals with respect to the input light 
modes flout ('^) = '^ini~^) with the complex conjugated 
spatial profiles. If the memory is operated off-resonantly, 
the optimal light modes and spin-waves may have differ- 
ent transverse profiles from the ones presented here, since 
both the absorption and refraction during light propaga- 
tion will be affected. 



late the beam waist and the distance between the focal 
plane and the beginning of the ensemble for the optimal 
input light modes. 

The results presented here are of direct importance for 
the experiments currently performed, where both the op- 
tical depth do and the Fresnel number F are limited by 
practical constraints. It may serve as a guide for choos- 
ing the optimal input light pulses to store given a set of 
physical parameters F and do , which may be constrained 
by, for instance, the limit in the achievable density (op- 
tical depth) of cold atoms or imperfection of the initial- 
ization by the optical pumping. Thus the knowledge of 
which light modes have optimal efficiency for such quan- 
tum memories is crucial [ISi] . 
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Appendix A: Proof for optimal read-out 



V. CONCLUSION 

We have solved the three-dimensional problem of light 
storage in a A-type atomic ensemble of finite spatial ex- 
tent. Specifically, we considered cylindrically symmetric 
samples with a constant atom density along the axial 
direction z. We have shown that within the paraxial ap- 
proximation, the only two important physical parameters 
characterizing the memory are the peak optical depth do 
and the Fresnel number F of the atomic ensemble. The 
maximal efficiencies for retrieval only as well as the full 
process of storage followed by retrieval have been calcu- 
lated. We considered both the forward and the backward 
operation modes of the memories and showed that both 
lead to high efficiencies for dense atomic media do > 40 
already for small Fresnel numbers F > 0.1. For small 
Fresnel numbers it is favorable to operate in the for- 
ward direction, whereas for large F, backward read-out 
yields higher efficiencies. The optimal incoming modes 
and corresponding stored spin-waves and output light 
modes have been presented for a set of physical parame- 
ters F = 2 and do = 100 for both backward and forward 
operating memory. The purity values of the optimal in- 
coming light modes are close to one, which is advanta- 
geous for experimental realizations. Moreover, we calcu- 



In this Appendix, we present the proof that the read- 
out efficiency is independent of the detuning A and the 
temporal shape of the control field 

We apply the spatial Laplace transform to the equa- 
tions of motion p6 | - (|28p and we find that the relation 
in Eq. (|34|) corresponds to the fact, that each matrix 
A{u',u) should obey the Sylvester matrix equation 




A{u' , u) 



(Al) 



where I is the unity matrix and the coupling matrix B, 
Eq. [20], has been written in the compact matrix nota- 
tion. We can write this equation as LA + AR — —I with 



L 



4irF 



■)- 



and 



Id, 



oK(u 



u '^^p ) B. It has a solution when the eigenvalues of the 
two matrices L and R obey: EIG(L) 7^ -EIG(R) [H, 
where EIG denotes the set of eigenvalues. We are only 
interested in purely imaginary u and u' , where u = iy 
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and u' = iy' , then Eq. (jAip reads 



Here, the matrices are 




A(m', u) 



(A2) 



+Hu',u)\\-^ 



The couphng matrix B is real and symmetric. Thus, one 
can see that the eigenvalues of the matrix L fulfill 



EIG(L) = EIG ( - 



y' 



= - +iEIG(W'), 



AttF 



where W = jd^)E> [y' — B is a real and sym- 

metric matrix. On the other hand the eigenvalues of the 
matrix R can be written as EIG(M) = i - iEIG(W) with 

r2 2 ■ ^ 



a real matrix 



Thus, the real 



parts are always equal, Re[EIG(L)j = Re[EIG(R)] = \. 
In consequence, the condition EIG(L) ^ — EIG(]R) al- 
ways holds, which proves that there is a solution to the 
Sylvester equation (|A1[) . This means that one can rewrite 
the position dependent loss and thus the read-out effi- 
ciency in terms of only the matrix A(m', u) and the spin 
wave 5'(tt, 0), which only depend on the optical depth do 
and the Fresnel number F while being independent of 
the detuning A and the temporal shape of the control 
field n{t). 



Appendix B: Alternative approach to forward 
read-out 



In Section IIV Al we presented the analytical expres- 
sions for the efficiencies of the storage followed by for- 
ward read-out, where Laplace transform in space was 
employed. In this Appendix, we present an alternative 
approach, where we perform Laplace transform in time, 
which for the storage procedure is defined by C{g[t)} = 
/o°° e~'^*g(t)(if, where w = w/7 is the dimensionless fre- 
quency. For the retrieval, the integral limits are from ~oo 
to 0. Here, for simplicity we assume the driving field to 
be constant in both space and time r2(r, t) = (see dis- 
cussion in Sees. HIl and IIV A]) and perform the derivation 
also in the co-moving frame. The input-output relations 
between the light modes and the stationary excitations 
then read 



flout (w) = / dzK[Vt,Cj,z\ Sq{z), 



(Bl) 



So{z) = — / d£oK'm*,Cj,l-~z\din{Co).{m) 



H(c^), 



E(w) = -i 



47rF 4J)(iA + i) 



If we normalize the incoming light mode according to 
2^1 /-ilo doj\din{u})\'^ = 1, the storage efficiency is given 
by 



1 

4^ 



dz I duj 

—i-co J — l-OO 



du'dljyCS) (B3) 
xJll"[r2*, w, 1 - z]K''[f7*, w, 1 - z]ain(w )■ 



The relation between the incoming and outcoming 
light modes is 



flout (w) = 7r~ / dz 
2m Jo 



dw'K[f],<I;,5] (B4) 



xK'^[j]*,w',l-z]ain(w') 

and thus, the combined efficiency for storage followed by 
forward retrieval is 



1 



Vs+h 



i-oo pi pi pi- 00 pi- 00 

dCj I dz I dz' I du I du' 

-i-oo ^0 ^ J — i-oo J —i-00 



xdl{iy)K*[n*,i^, 1 - ijKi'ifi,!:;), s] 



(B5) 



Analogously to the procedure discussed below Eqs. (jH]) 
and (|46|) . we diagonalize the kernel matrix ~ K*[f2*, ly,! — 
z]Kt[fi,c:;,z]]K[f^,t;),i']K^[fi*,z/',l - 5'] in order to find 
the eigenvalues giving the efficiencies of the storage fol- 
lowed by forward read-out. The corresponding eigen- 
vectors represent the optimized set of incoming light 
modes ain(0, v). 



Appendix C: Read-out in the backward direction 

In this Appendix, storage followed by backward re- 
trieval is analyzed. In the ID model, the optimal spin- 
wave for retrieval is real and the optimal storage effi- 
ciency follows from the optimal read-out efficiency by 
means of a time reversal argument as simply rj^ — For 
the three-dimensional problem, the optimal spin-wave is 
in general complex and in order to use the same argu- 
ment, one would have to physically perform complex con- 
jugation of the spin-wave, which is impossible. Thus, we 
need to optimize the full procedure as in the case of the 
forward memory. 

After we store the light mode in the spin- wave So{z) we 
want to retrieve it with a control field propagating in the 
opposite direction to the one used for storage. One way 
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to solve this problem is to rewrite the equations of mo- 
tion ([57| and for the opposite propagation direction 
(— z). On the other hand, one can also flip the spin- wave 
according to So{z) — >■ 5*0(1 — z) and next perform the 
read-out in the forward direction. In the Laplace trans- 
formed spatial variables, this transformation translates 
So{u) into 5o(-ii)e-™. 

Note that this method of describing backward read-out 
by the equations for forward propagation is only valid 
for degenerate ground states |0) and |1). Otherwise, 
the momentum conservation would be broken and the 
proper spin-wave to retrieve from in the backward direc- 
tion would read 5*0(1 — z)e^*'^'^^ with Ak = L{ujq — uji)/c. 
This accumulated phase leads to a decrease in the cffi- 
ciency of the memory operating in the backward direc- 
tion In Ref. [l^, it was proposed to partially over- 
come this problem by phase matching with signal and 
control pulses that are not collinear. 

In contrast to the forward memory, one needs to trans- 



form the spin-wave back to position space z after the 
storage is completed, otherwise the read-out is carried 
out not only from the atomic ensemble but also from the 
leaked excitation with z > 1. The input-output beam 
splitter relations for the backward read-out are 



aout(f) = ^ 



du / dze 

OQ Jo 



M[n{i),i,-u] 5o(z), 



So{z) 



1 



dt 



du e 



[n*i-i),-i,u] ai„(i). 



giving the outcoming light mode 



doutii) 



dt' 



du 



du' 



gi{u'-u) _ 2^ 

-ATi'^i{u' — u) 



xmSl{i),i,~u]M^p*{~i'),-i',u'] d;n{i'). 



Thus, the final expression for the efficiency of storage 
followed by backward retrieval reads 



f?s+br = 



oo poo 



— oo 'J — oo <J —oo 'J —oo <J —oo 'J —oo 



oo poo 



di^ I dt I di' I du I du' I du I du'd\^{i) 



^-i{u'-u) _ ^ 

4:'K'^i{u' — u) 



X M*[n*{-i),-i,u'] M''[n{i^),i^,-u]M[n{Q,i^,-u] M^'[n* {-i'),-i',u 



(CI) 



-4Tr^i{u' — u) 



am(?) 



In the position space and frequency picture, the re- 
lation between the incoming and outcoming light mode 
is 



dCb' Kin^Cj^S] (C2) 



aout(w) = — / dz 



xIK^[17*,c;)',z]ain(w'), 

which leads to the final expression for the memory work- 
ing in the backward direction 



1 



l-OO pl-OO 



r-l rl 

duj I dz I dz' I diy I du' 



^ J -i-co Jo Jo J—i-oo J —i-oo 

xdl{i^)K*[n*,i^,s]K''[n,Q,s] (C3) 

xK[0, u, S']K^[n*, v', z']d;^{v'). 

Again, we optimize the memory by looking at the 
eigenvalues and eigenvectors of the kernel matrix 
K*[h\v,i]K'^[n,Cj,2]K[n,Co,i'\K^[n\v',i'] or 
- M* -i, u']Mt [fi(i), i„ -w]M[f2(i), L, -ii] 

M.^[Cl* {—i'), —i' ,u'], which gives us the efficiencies for 
the storage with backward read-out and the optimized 
incoming light modes, respectively. 



Appendix D: Proof of the relation between the 
input and output light mode 

In this Appendix, we prove that the output light is 
the time reversed complex conjugate of the input mode 
flout (0 — 0*11 (~^) fo'^ both backward and forward oper- 
ating memories. We start with the forward memory and 
the beam splitter relation between the input and output 
light modes in the adiabatic limit with a constant Rabi 
frequency and make a substitution of the integral time 
variable t' ^ —t' leading to: 



flout (i) 



du 



di' — 

2TTi 



f^[i',u] ai„(-?). 



at the transpose 



By looking 

trix Bi[i,i'] J^^duM[i,u]M 



of the ma- 

M^[i', 



, u 



see that it 



is symmetric Df[i,t'] = Df [t',i\. Singular value 
decomposition of this matrix can be written as 
D[i,t'] = \/Vi{''Ji\ with the square roots of the 

corresponding efficiencies. The fact that this matrix is 
symmetric implies that \ui) — \v*), which means that the 
output light is the time reversed complex conjugate of 
the input light mode, flout (0 = '''rn(~*)- transverse 
intensity shape of the output mode is the same and the 
phase is the complex conjugate of the input light mode. 
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Similarly for the backward operating memory, one can 



substitute t' ~t' and 



and write the beam- 



splitter relation in Eq. (jClj) as 



flout (i) 



di' 



du 



du' 



-47r2i(u' + u) 



X 



\u]M^[i',u'] ai„(-f')- 



Again, we consider the transpose of the ma- 
trix Db[t,?] = J^^duJ^^du'U[i,u\M^[i',u'] ^ 

/^dw/^dw'Q(w)ef*(")*eN''("')*'Q'r(M'), which reads 

Substituting u ^ u' , we see that this matrix is sym- 
metric, Db[t, ?] = Db^[i',f] so that also in the case of 
backward read-out we have the above symmetry between 
the input and output light modes. 



Appendix E: Influence of the purity on the 
efficiencies 



As we have shown in Sec. lIVBl the optimal input light 
modes are not necessarily pure, however, the impurity 
is small and does not go above 4.5% for the parameters 
considered here. Experimentally, one would prefer to use 
pure light pulses to avoid the difficulty of changing the 
transverse profile in time. Similarly, it is preferable to 
detect only the dominant Schmidt component of the cor- 
responding output light. In the main text, we present effi- 
ciencies of the fully optimized pulses without constraints 
on light pulse purities. In the present appendix, we es- 
timate an upper bound to how much the experimental 
simplifications of using separable modes hn{r±) and a„(t) 
will affect the efficiencies. 

One can write the optimal input and output light 
modes as a sum of its dominant Schmidt component 
^hi/out (^) = "fn*/outffin/out(^) and the remaining impure 
part (correlated in time t and transverse posi- 

tion f±) 



Tit) = /f^a P(t) + v/T^^^(i), 



(El) 



* opt 



respectively. In the space of input light modes, we now 
consider the 2-dimensional subspace Lin spanned by the 
ordered orthonormal basis /3in = {«fn(^)i fl|n(^)}- Ob- 
viously a°f*(?) G Lin with the coordinate representa- 
tion ei = ( /T—^j-^- We then define d^{i) £ Lin, 

as the mode with co- 



which is orthogonal to a°n'(?) 
ordinate representation 6*2 



Pi, 



Now we 



consider the map G describing storage followed by re- 
trieval restricted to the domain Lin (G is an abbrevia- 
tion for input-output relations of the type that appear in 
Eq. (HSI) ): denote the range of G by Lout- By definition. 



a°f*(i) is an eigenmode of G'^G with eigenvalue ?7max and 
we have ^?7maxa^ut(^) = '^«in*(^)- We define the nor- 
malized mode d^^^{t) by Vrf'^outii) — ^^'^ since 
«hi(*~) -L a°f'(t), we have 

«ou\(i) • cut) = -^==[aiHt)G^] ■ [Ga°r(t)] 

V '7max'7 



V '/max'/ 

= V " 

that is, d^^^{i) _L ^^^(f); here • signifies the inner prod- 
uct in Bessel mode space and time. Next, we intro- 
duce an ordered orthonormal set /3out of output light 
modes, /3out = {d'^^^{i),d'^^^{t),b{i)}, where b{t) is a 



normalized mode that fulfills cb{t) = d^^^.{t) — [a^i^tit) ■ 

a^ut(?)]aout(t) - [aout(t) ' aiit(t)]aout(t~), for some c G C. 
We have Lout ^ span(/3out) and the coordinate repre- 
sentations w.r.t. /?out of a^ut(i) and a^ut,(?) are /i = 

V-Pout \ ^ / ~V1 ~ Po' 

and /2 = 



\/l ^ -Pout 



spectively, for some 



V Pout 

c' 



re- 



G M, c' G C; the ratio of the first 
two coordinates of /2 follows from a^^^{i) _L aTt{i). 

We can now deduce a lower bound for the efficiency 
'7puro of storing the pure mode af^ (t) followed by retrieval 
where the output light is projected into the spatial pro- 
file a^ut of the pure mode^aPnt(t) = aout3out(t~)- Note, 
however, that in general a^^^ is not the spatial profile 
with the largest integrated overlap with the output light 
corresponding to df^(t). In the coordinate space of Lin, 

df^{i) is represented by ( ^ ) = ^/^el - \/r~~^e2. 



whereby 



77purc = 1(10 0) [G]^™' I 

= 1(1 ) (V^/f^/i - VrfVl - Pinf2 



I \/ '7max \/ Pin \/ Po\ 



?/Vi - p'hVi - Po, 



Vi + FF 



Assuming that the largest eigenvalue Ty^ax of G^G is non- 
degenerate (which seems to be generally true from nu- 
merical calculations [31]), we have from d^{i) _L a°n*(t) 
that rj' < ri2, where 772 is the second largest eigenvalue. 
However, since only the values of f^max are presented in 
the present paper, we will assume rj' — jymax resulting in 
a smaller lower bound. Maximal destructive interference 
will occur for (f) = tt and c' = 0; making these substitu- 
tions along with Pout = Pn (as follows from the result of 
Appendix ID|) . we obtain 

?ypurc > '7max(l " 2Pi„)^ = ?7max(l " 2e)^ W ?ymax(l ^ 4e), 
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where e = 1 — Pin is the impurity. The efRciency of the 
full optimal mode in the worst case of P = 0.9581 for 



forward retrieval, to = 0, do = 200, and F = 0.02 can 
at most be reduced by ~ 16.76% (from r] = 0.8049 to 
?7pure = 0.67). 
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